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Abstract 

The values of ion-atom ionization cross sections are frequently needed for many applications 
that utilize the propagation of fast ions through matter. When experimental data and theoret- 
ical calculations are not available, approximate formulas are frequently used. This paper briefly 
summarizes the most important theoretical results and approaches to cross section calculations in 
order to place the discussion in historical perspective and offer a concise introduction to the topic. 
Based on experimental data and theoretical predictions, a new fit for ionization cross sections is 
proposed. The range of validity and accuracy of several frequently used approximations (classical 
trajectory, the Born approximation, and so forth) are discussed using, as examples, the ionization 
cross sections of hydrogen and helium atoms by various fully stripped ions. 



1 



I. INTRODUCTION 



Ion-atom ionizing collisions play an important role in many applications such as heavy 



ion inertial fusion l|, coUisional and radiative processes in the Earth's upper atmosphere 

n fi 

ion-beam lifetimes in accelerators 3], atomic spectroscopy and ion stopping in matter 
S. and „e of —able M.e. a.o.,c p.... g T.e .ce. resurgence of interest in 
charged particle beam transport in background plasma is brought about by the recognition 
that plasma can be used as a magnetic lens. Applications of the plasma lens ranging from 
heavy ion fusion to high energy lepton colliders are discussed in Refs. [6-10]. In particular, 
both heavy ion fusion and high energy physics applications involve the transport of positive 
charges in plasma: partiall;j;_^ stripped heavy elements for heavy ion fusion; positrons for 
electron-positrons colliders and high- density laser-produced proton beams for the fast 
ignition of inertial confinement fusion targets 11 1. 

To estimate the ionization and stripping rates of fast ions propagating through gas or 
plasma, the values of ion-atom ionization cross sections are necessary. In contrast to the 
electron 2| and proton , Q] ionization cross sections, where experimental data or 

theoretical calculations exist for practically any ion and atom, the knowled ge o f ionization 
cross sections by fast complex ions and atoms is far from complete 

0,0, HQ- When 

experimental data and theoretical calculations are not available, approximate formulas are 
frequently used. 

The raison d'etre for this paper are the frequent requests that we have had from colleagues 
for a paper describing the regions of validity of different approximations and scaling laws 
in the calculation of ion- atom stripping cross sections. The experimental data on stripping 
cross sections at low projectile energy were collected in the late 1980s, while comprehensive 
quantum mechanical simulations were performed in the late 1990s. Having in hand both 
new experimental data and simulation results enabled us to identify regions of validity of 
different approximations and propose a new scaling law, which is the subject of the present 
paper. 

The most popular formula for ionization cross sections was proposed by Gryzinski {2^. 
The web of science search engine j2l| shows 457 citations of the paper, and most of the 
citing papers use Gryzinski's formula to evaluate the cross sections. In this approach, the 
cross section is specified by multiplication of a scaling factor and the unique function of the 
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projectile velocity normalized to the orbital electron velocity. The popularity of Gryzinski's 
formula is based on the simplicity of the calculation, notwithstanding the fact that his 
formula is not accurate at small energies. 

Another fit, proposed by Gillespie, gives results close to Gryzinski's formula at large 
energies, and makes corrections to Gryzinski's formula at small energies |22]. Although 
more accurate, Gillespie's fit is not frequently used in applications, because it requires a 
knowledge of fitting parameters not always known a priori. 

In this paper, we propose a new fit formula for ionization cross section which has no fitting 
parameters. The formula is checked against available experimental data and theoretical 
predictions. Note that previous scaling laws either used fitting parameters or actually did 
not match experiments for a wide range of projectile velocities. We also briefiy review the 
most important theoretical results and approaches to cross section calculations in order to 
place the discussion in historical perspective and offer nonspecialists a concise introduction 
to the topic. 

The organization of this paper is as follows. In Sec. II we give a brief overview of key 
theoretical results and experimental data. Further details of the theoretical models are 
presented in Appendices A-C. The new proposed fit formula for ionization cross section is 
presented in Sec. Ill, including a detailed comparison with experimental data, and in Sec. IV 
the theoretical justification for the new fit formula is discussed. 



II. BRIEF OVERVIEW OF THE THEORETICAL MODELS AND EXPERIMEN- 
TAL DATA FOR IONIZATION CROSS SECTIONS 

There are several theoretical approaches to cross section calculations. These include: 
classical calculations that make use of a classical trajectory and the atomic electron velocity 
distribution functions given by quantum mechanics [this approach is frequently referred to 



as classical trajectory Monte Carlo (CTMC)]; quantum mechanical 



calculations based on 



Born, eikonal or quasiclassical approximations, and so forth All approaches 

are computationally intensive and the error and range of validity are difficult to estimate in 
most cases. Therefore, different fittings and scalings for cross sections are frequently used 
in practical applications. 

Most scalings were developed using theories and simulations based on classical mechanics. 
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Classical trajectory calculations are easier to perform compared with quantum mechanical 
calculations. Moreover, in some cases the CTMC calculations yield results very close to 
the quantum- mechanical calculations 0, 2h) 22 1. The reason for similar results lies 



in the fact that the Rutherford scattering cross section is identical in both classical and 
,uant.„-.echa„,cal denvat.o,. Q. The^efo.. whe. an ion,..g colHs.on . predominantly 

a consequence of the electron scattering at small impact parameters close to the nucleus, 
the quantum mechanical uncertainty in the scattering angle is small compared with the 



angle itself, and the classical calculation can yield an accurate description 
Whereas in the opposite limit, when an ionizing collision is predominantly a consequence 
of the electron scattering at large impact parameters far from the nucleus, the quantum 
mechanical uncertainty in the scattering angle is large compared with the angle its elf, and 
the classical calculation can remarkably fail in computing the ionization cross section (3ll.l32j|. 

In the present analysis, we consider first the stripping or ionization cross section of the 
hydrogen- like electron orbits (for example one-electron ions), with nucleus of charge Zt, 
colliding with a fully stripped ion of charge Zp. Subsequently, we show that the approach 
can be generalized with reasonable accuracy for any electron orbital, making use of the 
ionization potential of the electron orbitals. Because different terminology is used in the 
literature, we call a stripping collision a collision in which the fast ion loses an electron in 
a collision with a stationary target ion or atom (in the laboratory frame); and we call an 
ionizing collision a collision in which a fast ion ionizes a stationary target ion or atom 3| ■ 
Both cases are physically equivalent to each other by changing the frame of reference, and 
further consideration can be given in the frame of the atom or ion being ionized. 

Atomic units are used throughout this paper with e = h = nif. = 1, which corresponds 
to length normalized to = ^^/(meC^) = 0.529 ■ 10~^cm, velocity normalized to vq = 
e'^/h = 2.19 ■ lO^cm/s, energy normalized to Eq = nieVQ = 2Ry = 27.2eV, where Ry 
is the Rydberg energy. The normalizing coefficients are kept in all equations for robust 
application of the formulas. For efficient manipulation of the formulas it is worth noting 



that the normalized projectile ion velocity isv/vo = Q.2y E[keV / amu], where is energy per 
nucleon in keV/amu. Therefore, 25keV/amu corresponds to the atomic velocity scale. Some 
papers express the normalized velocity v/vq as Pa, where /? = v/c, and Vq/c = a = 1/137. 
Here, c is the speed of light, and a is the fine structure constant. 

For a one-electron ion, the typical scale for the electron orbital velocity is Vni = VqZt- 
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Here, n, / is the standard notation for the main quantum number and the orbital angular 
momentum The collision dynamics is very different depending on whether v is smaller 
or larger than v^i- 



A. Behavior of cross sections at large projectile velocities v > Vni 

If V >> Vni, the electron interaction with the projectile ion occurs for a very short time 
and the interaction time decreases as the velocity increases. Therefore, the ionization cross 
section also decreases as the velocity increases. In the opposite case v « Vni, the electron 
circulation around the target nucleus is much faster than the interaction time, and the 
momentum transfer from the projectile ion to the electron averages out due to the fast 
circulation. Thus, the cross section decreases as the projectile velocity decreases. This is 
why the cross section typically has a maximum at f = fmax ~ Vni, but as we shall see below, 
Vmax also dcpcuds on the charge of the projectile. 



1. Thompson's treatment 

In the first treatment, Thompson calculated the ionization cross section in the limit 
V >> Vni 01 . This treatment neglected completely the orbital motion of the target electrons 
and assumed a straight-line trajectory of the projectile. In this approximation, the velocity 
kick acquired by the electron during the collision is entirely in the direction perpendicular 
to the ion trajectory, because the final action of the force along the trajectory cancels out 
due to symmetry, i.e., the electron velocity change during the approaching phase is equal 
to minus the electron velocity change during the departing phase. The momentum acquired 
by the electron ( rrieAv) from passing by the projectile moving with the speed v and impact 
parameter p is given by the integral over time of the force perpendicular to ion trajectory 
F± = e^Zpp/ (p^ + f ^t^)^/^, where t = corresponds to the distance of the closest approach. 
Time integration of the force yields 

Avip) = (1) 
rrieVp 

From Eq. (Q) it follows that only collisions with sufficiently small impact parameters result in 
ionization. The minimum impact parameter for ionization of an initially stationary electron 
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(Pmin) is mgAf (pmin)^/2 = Duimg a collision with impact parameter pmin the energy 
transfer from the projectile to the electron is equal to the ionization potential /„/ = Z^Eq/2, 



or Av 



33| 



Pmin) = Vni- Substitution of Eq-iH)) gives the total ionization cross section vrp 



2 

min 



^Bohr I 



2 2 v'^Eq 



U 1; 



nl 



(2) 



Similarly, Eq.Q can be derived by averaging the Rutherford cross section over all scattering 
ang les leading to ionization. Although the first derivation of Eq.( 21) was done by Thompson 
33[ the formula is frequently referred to as the Bohr formula 



2. Gerjuoy's treatment 

The following treatments account for the effect of finite electron orbital velocity. The 
most complete and accurate calculations were done by Gerjuoy 3^. He calculated the 
differential cross section da/dAE{ve,v, AE) of energy transfer AE in the collision between 
the projectile ion and a free electron (the target atomic potential was neglected) with given 
initial speed Ve (and arbitrary direction), by averaging the Rutherford cross section over 
all orientations of electron orbital velocity Vg. The total cross section is then calculated by 
integration over the energy transfer for energies larger than the ionization potential, and 
weighted by the electron velocity distribution function / (t>e). This gives 

POO 

a{v,InhZp) = Zp ai^^X^.Ve) f {Ve) dVe, (3) 

where 

-^iv,Ve,AE)dAE. (4) 

A rather complicated analytical expression for da /dAE{ve,v, AE) is given in Appendix A. 
For large projectile ion velocities {v >> Vni), the differential cross section can be expressed 
as 3 

, high-energy /r, 2 \ 

Substituting Eq.® into Eq.Q and Eq.(jll) gives 

"S.ZT''^. h„ Z,) = jB„,^'"""-(i,, /,.„ Z,), (6) 



nl 
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where a^°^^is given by Eq.Q, and Kni =< >ni is the average orbital electron kinetic 

energy. For hydrogen-like electron orbitals, the average electron kinetic energy is equal to the 
ionization potential Kni = Ini [s^]' ^ni = 1- The Bni factors are introduced to account 
for the difference in the electron velocity distribution functions (EVDF) from the EVDF of 

ihe hydrogen-like electron orbitals. The data for K^i are calculated for many atoms in Ref . 

3^. For example, the average kinetic energy for the helium atom is Kni =< >= 
1A3Eq, whereas /„; = 0.91-Eo; and therefore Bhb = 1-22. That is the reason that accounting 
for the finite orbital electron velocity gives a cross section which is 5/3 times larger than the 
Bohr formula in Eq. . This is a consequence of the fact that for an electron with nonzero 
velocity less energy transfer is required for ionization. 

Classical mechanics gives the EVDF as a microcanonical ensemble, where 



Here, C is a normalization constant defined so that J / (tig) dVf. = 1, and denotes the 
Dirac delta-function. Interestingly, the EVDF for a one-electron ion is identical in both the 
,uantu™chanic^ and c.a..a, calculation 0Qw.H 

where Vni is the scale of electron orbital velocity 



Vnl = f 



V^hjE~o- (9) 



Although a microcanonical distribution provides the same velocity distribution as in quan- 
tum theory for hydrogen-like shells, this is not the case for other electron shells. Moreover, 
the spatial distribution of the charge density is poorly approximated even for hydrogen, 
vanishing identically for r > 2ao rather than decreasing exponentially . Substituting the 
general differential cross section da/dAE{ve,v, AE) from Eq. ()A3|) of Appendix A and the 
EVDF in Eq.® into Eq.© yields 

a^^^'iv, Z,) = T^alZl ^G^^"" (^) . (10) 

^nl \VnlJ 

Here, the scaling function G'~'^^{x) is given by Eq. ()C3|) in Appendix A, using the tabulation 



of the function G{x) presented in Ref.j36| for x > 1, and in Ref. 



GGV stands for the classical trajectory calculation in Eq. ()C3|) due to Gerjuoy 
fit of Garcia and Vriens 



for X <1. The notation 
using the 



3. Bethe's treatment 

The classical calculations underestimate the cross sections for very high projectile ve- 
locities V >> Vni- The scattering angle of the projectile due to collision with the target 
atom is of order 6^ = Ap/Mv, where Ap is the momentum transfer in the collision, and M 
is the mass of the projectile particle. The minimum energy transfer from the projectile is 
determined by the ionization potential, with AE = vAp > Inu and Ap > Apmin = Ini/v. 
Here, we use the fact that the momentum transfer Ap is predominantly in the direction 
perpendicular to the projectile velocity. The projectile particle with wave vector k = Mv/h 
undergoes diffraction on the object of the target atomic size ani with the diffraction angle 
of order 6^ = l/{kani) = h/{Mvani) large projectile velocities v >> Vni, it follows 

that Apmin = Ini/v « k/anh bccausc Vni = Inio.nl foT hydrogcu-like electron orbitals. 
And for small Ap ~ Ap^[^, it follows that 6c = Ap/Mv << 9^ = h/ (Mvani)- Therefore, the 
collision can not be described by classical mechanics. 

Bethe made use of the Born approximation of quantum mechanics to calculate cross 
sections js^ (see Appendix B for details). This yields for i; >> Vni 



If the projectile speed is much larger than the electron orbital velocity v >> Vnu the loga- 
rithmic term on the right-hand side of Eq. ^\\\ contributes substantially to the cross section, 
and as a result the quantum mechanical calculation in Eq. (fTT|) gives a larger cross section 
than the classical trajectory treatment in Eq.®. The quantum mechanical cross section 
is larger than the classical trajectory cross section due to the contribution of large impact 
parameters (p) to the quantum-mechanical cross section, where the ionization is forbidden in 
classical mechanics because the energy transfer calculated by classical mechanics is less than 
the ionization potential [AE = vApc{p) < Ini, where Apc is the momentum transfer given by 
classical mechanics in Eq.(^]. However, ionization is possible due to diffraction in quantum 
mechanics js^. Moreover, integration over these large impact parameters where the ion- 
ization is forbidden in classical mechanics, contributes considerably to the total ionization 
cross section (see Appendix B for further details). 




(11) 
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4- Gryzinski's treatment 



Gryzinski attempted to obtain the ionization cross sections using only classical mechanics 
similarly to Gerjuoy. But, in order to match the asymptotic behavior of the Bethe formula 
in Eq. ljllj) at large projectile velocities, Gryzinski assumed an artificial electron velocity 
distribution function (EVDF) instead of the correct EVDF in Eq.® [2^, i.e.. 



/^^^M^e) =-(^Vexp(-^). (12) 

Vnl \Ve J \ Ve J 

The ionization cross section was calculated by averaging the Rutherford cross section over 
all possible electron velocities, similar to the Gerjuoy calculation in Eq.®, but was less 
accurate for small velocities v < Vni- The effect of using the EVDF in Eq. p2|) is to 
populate the EVDF tail with a much larger fraction of high-energy electrons with Ve » Vni, 
which gives f^^^^ (ve) ~ instead of / (fe) ~ for the correct EVDF in Eq.(|H)l. 
As a result, the average electron kinetic energy < mf.vl/2 > diverges, which leads to a 
considerable enhancement of the ionization cross section at high projectile velocities. For 
V >> Vnl, Gerjuoy's calculation of the differential cross section da/dAE{ve, v, AE) of energy 
transfer AE is similar to Gryzinski's. Therefore, we can substitute Eq. fll2j) into Eqs.(j3|) 
and (j3)). Because in the limit v » Vni the ionization cross section is proportional to 
the average electron kinetic energy < mevl/2 > [Eq.^], and the average kinetic energy 
diverges, it follows that a small population of high-speed electrons contributes considerably 
to the cross section. Using the general expression for d(T/dAE{ve,v, AE) avoids singularity 
and yields the logarithmic term in the ionization cross section similar to the Bethe formula in 
Eq.(fTT|). After a number of additional simplifications and assumptions, Gryzinski suggested 
an approximation for the cross section in the form given by Eq. pU|) with 20| 

a^'^y^iv, hu Z,) = T^alZl ^G^'^y^ ( —) . (13) 

^nl \VnlJ 

Here, the function G'^^^^(x) is specified by Eq. flG6|l of Appendix C. In Eq. ()13|) . the function 
QGryz^^^ has the following limit 

G^'-y'ix) ^[l + 0M7ln{2.7 + x)]/x^ as x ^ 00, (14) 
which is close to Bethe's result in Eq. (fTT|) . 

QBethe^^-j [1.261 + 0.566 ln(a;)] /x^ as x oo. (15) 



For 10 < a; < 40, it follows that 



G^'y'{x)/G^^"'\x)^lM. (16) 

Therefore, the Gryzinski formula can be viewed as a fit to the Bethe formula at large 
velocities f >> Vni with some rather arbitrary continuation to small velocities v « Vni- 

Figure shows the experimental data for the cross section for ionizing collisions of fully 
stripped ions colliding with a hydrogen atom, 

+ H{ls) + H+ + e, (17) 

where X'^'^ denotes fully stripped ions of H, He, Li, C atoms, and (Is) symbolizes t he g round 
state of a hydrogen atom. The experimental data for ions were taken from 40| (note 
that authors of this reference concluded that the previous measurements of the cross sections 
were inaccurate); from for He^^, C^^ ions ; and from j3l for Li'^^ ions. 

From FigH it is evident that the Bethe formula describes well the cross sections for 
projectile velocities larger than the orbital velocities v » Vni- At large energies, the GGV 
formula underestimates the cross sections as discussed before, whereas Gryzinski's formula 
gives results close to the Bethe formula and the experimental data. Both, the GGV and 
Gryzinski formulas disagree with the experimental data at small energies. 

B. Behavior of cross sections at small projectile velocities v < Vni 

The Bethe, GGV and Gryzinski's formulas fail at small velocities because they assume 
free electrons, neglecting the influence of the target atom potential on the electron motion 
during the collision. Apparently the assumption of free electron motion fails if the circulation 
period of the electron around the atom's nucleus is comparable with the interaction time 
of an ion with the electron. Let us now estimate the projectile velocities at which the 
electron circulation needs to be taken into account. The typical impact parameter leading 
to ionization is 



TC VVnl 

and the interaction time is of order Pioniz/v. The electron circulation time is r„; ~ ani/vni, 
where Vni is the electron orbital velocity, which scales as Vni = ZtVq, and ani is the ion radius 
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FIG. 1: Ionization cross sections of hydrogen by fully stripped ions showing both experimental data 
and theoretical fits. GGV stands for the classical calculation in Ea. (|G3|) due to Gerjuoy using the 
fit of Garcia and Vriens. Gryz. denotes the Gryzinski approximation in Ea. ()C6() . Bethe stands for 
Bethe's quantum-mechanical calculation in the Born approximation, limited to u > Vni in Ea. Hllj) . 
Finally, BA denotes the Born approximation in the general case in Eq. (|2H) . All values are in 
atomic units. All values are in atomic units. For hydrogen, the ionization potential is = \/2Eq, 
Vni = vq = 2.19-lO^cm/s, and the cross section is normalized to tto^/I'^i = 3.51- 10~^^cm^. Symbols 
show experimental data. 

o-ni = (Iq/Zt IsQ^. Therefore the condition r„/ > pioniz/v holds for v > fmax, where 

Wmax = Vnl\j2Zp/ZT. (19) 

Here, Zp is the charge of the fully stripped projectile and Zt is the nuclear charge of the 
target atom or ion. For velocities larger than fmax, the ionization cross section decreases as 
the velocity increases [see Eq.llllj)] due to the decreasing interaction time with an increase 
in velocity. On the other hand, for velocities less than fmax? the collision becomes more 
adiabatic. The influence of the projectile is averaged out due to the slower motion of 
the projectile compared with the electron orbital velocity, and the ionization cross section 
decreases with decreasing projectile velocity. Thus, the cross section has a maximum at 

11 



max 



[Eq.dlHD]. 



Note that if the projectile speed is comparable with or smaller than the electron orbital 
velocity v < Vni, the Born approximation of quantum mechanical theory is not valid. Cum- 
bersome quantum mechanical simulations are necessary for an exact calculation of the cross 
sections, as for example in Ref. 43[. Nevertheless for the case 2Zp ~ Zt the maximum 
of the cross section calculated from the Born approximation is similar to the experimental 
results. To describe the behavior of the cross section near the maximum, the second-order 
co„.ect,o. ,n the p«te. been calculated ,n Q. yie.d.g the c.oss sect.o. ,n 

the form 

a^'l'^'iv) = ^ ^ 0.566 In f?^) + 1.26 - 0.66^ , (20) 



0.566 In (w) + 1.26 - 0.66^ 



y2 



'mod K^J j2 ~2 

where v = v/vni- Equation ()20|l agrees with the exact calculation in the Born approximation 
[Eq. ()Bl|) ] as described in Appendix B (the agreement is within 10% for t; > 1.1). We have 
developed the following fit for the cross section in the Born approximation in the general 
case, 



2 Z^ 



1.95 



v(l + 1.2v^) 



(21) 



af^fiv) = 9^ 4 [0.283 In {v' + l) + 1.26] exp 

Equation ()21|) agrees with the exact calculation [Eq. ()Bl|) ] within 2% for -y > 1, and within 
20% for 0.2 <v< 1. 

Equation (j2T| was derived making use of the unperturbed atomic electron wave functions, 
which implicitly assumes that the projectile particle transfers momentum to the electron and 
departs to large distances, where it does not affect the electron to be ionized. The wave 
function can therefore be described as a continuous spectrum of the atomic electron, not 
affected by the projectile. 

This assumption breaks down at low projectile velocities when the projectile velocity is 
comparable with the electron orbital velocity. Indeed, the electron kinetic energy in the 
frame of the projectile is of order mef^/2 and the potential energy Zpc'^ / pioniz, where pioniz 
is the impact parameter leading to ionization, given by Eq. (jl8p . Substituting pioniz from 
Eq. lfTHj) into electron potential energy Zpc'^/pioniz gives that potential energy is larger than 
kinetic energy if 

V<Vnl. (22) 



Therefore, under the condition in Eq. ()22jl . an electron can be effectively captured by the 
projectile after the collision instead of leading to ionization. As a result, the ionization cross 
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section is small compared with the charge exchange cross section at low projectile velocities. 
The assumption of the unperturbed electron wave function results in grossly overestimated 
ionization cross sections as can be seen in Fig.l. 

The ionization cross sections are also difficult to measure at small projectile energies, 
because careful separation between the large charge exchange cross section and the small 
ionization cross section is necessary for the correct measurement 0]. Therefore, early 
measurements of the ionization cross section at small velocities were not always accurate 



1. Gillespie's treatment 

To account for the difference between the Born approximation results and the experimen- 
ial data for v < fmax, Gillespie proposed to fit the cross sections to the following function 



Here, A„/ is a constant, which characterizes the ionized atom or ion (for example, for the 
ground state of H, Xni = 0.76), and cr^f^'^ is the cross section in the Born approximation in 
the form of Eq.()20j). Gillespie's Eq.()23|) proved to fit very well existing experimental cross 
sections for hydrogen atom ionization by He^^, Li^'^,Li~^^, C"*"^, N^^, N~^^, O^^ ions, 
and less well for He and H molecules with the same ions [s^- Because (J^fJ^^i^v) becomes 
negative for v < 0.7, Gillespie's Eq.(|^ can not be applied to these low projectile velocities. 
In principle, the general fit af^f in Eq. ipTj) can be used instead of c^f^'^ in Eq. (p]j) . However, 
because the two formulas differ considerably in the range of interest, 0.7 < f < 1, the fitting 
coefficients have to be updated for use with cr^-f. 

Although Gillespie's fit proved to be very useful, there are a number of reasons to look 
for another fit. Gryzinski's Eq. ()(y6|l is frequently used, because it requires only knowledge of 
one function for calculations of cross sections, notwithstanding the fact that it overestimates 
the cross sections at low energies. 



lfl|40|. 
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2. Bohr and Linhard's treatment 



For V ^ Vni, a. universal curve is expected if both the cross sections and the square of 

r~i 

impact velocity are divided by Zp [48]. This scaling was established for the total electron 
loss cross section o"*^', which includes both the charge exchange cross section a'^'^ and the 
ionization cross section. Based on the results of classical trajectory Monte Carlo (CTMC) 



calculations, Olson proposed the following fit 



49| 



a^\v, Zp) = ZpA^iT^alf'^"- — , (24) 



where /(x) describes the scaled cross sections 

Here, 7„; and Ani are constants, for example, '~)h = a/5/4 = 1.12 and Ah = 16/3 for 
atomic hydrogen, and = 1-44 and A^e = 3.57 for helium. The scaling in Eq. ()24|l 
was also demonstrated analytically by Janev j^]. For v « v^^fZp, a^^ is dominated by 
charge exchange, o""^^ ~ cr^', and Eq. ()24|l gives a constant cross section for charge exchange, 
cr'^*^ a^^ = WnZp/SaQ. For v >> v^^fzFp^ a^' is dominated by the ionization cross section, 
and a'^'^ ~ ^^dassicai'^^^ [Eq-®]- Note that the scaling in Eq. lj^ does not reproduce the 
logarithmic term in the Bethe formula [Eq. ljlll) ] for v » v^^^fZp because it is based on 
classical trajectory calculations. To make Eq. (f2l|l agree with Eq.®, the coefficients 7^/ 
should be proportional to \fT^i. For example, the ionization potential for hydrogen is Ih = 
13.6el^, and for helium Ihb = 24.6eV. The ratio of 7^ = 1.12 to 7^6 = 1-44 differs 
from y/Tfj/y/lHe by only five percent, i.e., ^nl^fhil {iHe/ V I He) = 1-05. Therefore, as was 
shown by Janev 1501 , the scaling in Eq. can be rewritten in a form similar to Eq. pop by 
normalizing the velocity to v^u Eq-©, i-e., 

a'\v, Ini, Zp) = T^alZpN^i ^Br^iG^' I ] , (25) 

^nl \Vnl\/Zpj 

where 

G^' (x) = ^/^^^°" {xhu) . 

Here, Nni is the number of electrons in the orbital n/, and the B^i factors Eq.© are in- 
troduced to account for the difference of the orbital electron velocity distribution functions 
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with the hydrogen-hke EVDF function in Eq.(jH)). By construction, Eq.()25p coincides with 
Eq.® in the hmit v » Vni\fZ'p- 

Because the scahng in Eq. lj^ is based on classical trajectory calculations, it is valid 
only for intermediate velocities where the underbarrier transitions allowed in the quantum 
mechanical calculations do not contribute significantly (see Appendix B for details). Experi- 
mental data 42, 5^ confirm the scaling in Eq.p^ for 1.2 < v / {vni\f^) < 3, or equivalently, 



for the projectile energy in the range = 30 — 200 x Zplni/In in units of keV/amu. 



A similar scaling to Eq. ()24|) was derived in Ref.jGJ] based on quantum mechanical cal- 
culations making use of the quasi-classical approach developed originally by Keldysh for 
multi-photon ionization of atoms in a strong electromagnetic field. These calculations give 
scaling similar to Eq. ()24j) . but with a different function f{x) given in The quantum 

mechanical calculation results for the charge exchange cross section in Ref. Q] are a factor 
of 3 larger than Olson's cross section in Eq. ()24j) for v j [y^^^fZ'^ < 0.2. 

Direct application of the scaling in Eq. ()25|l for the ionization cross section instead of 
the total electron removal cross section does not produce a single scaled function [see Fig. 2 
for hydrogen and Fig. 4.(b) for helium]. Furthermore, the data are considerably scattered 
near the maxima of the cross sections. A number of other semi-empirical models have been 
developed, which use up to ten fitting parameters to describe the ionization cross sections 
over the entire projectile energy range 19^ . 

III. NEW FIT FORMULA FOR THE IONIZATION CROSS SECTION 

Analysis of the experimental data in Fig^ shows that the maxima of the experimentally 



measured cross sections occur at Zp + 1, not at ^/ Zp as would be the case according to 
Olson's scaling in Eq. ()24|l . Therefore, it is natural to plot cross sections as a function of 
the normahzed velocity v/ {vni^/Zp + 1). Note that at large velocities, according to Eq.ijH} 



a ~ Zp/v . Therefore, making use of the normalized velocity v/ {vni^/Zp + 1) requires 
normalization of the cross sections according to a/ [Zp/{Zp + 1)] . As a consequence, instead 
of Eq. (|^^ , we propose the following scaling 

a-{v, Zp) = <n^^.N^i^G-^- [ ^==\ . (26) 



15 



Exp. □ o He V Li'^ a c 

1 



+6 



N 



^ 0.1 



1 1 1 1 1 1 

WW 
O V 






A E„ 

□ n » 

□ 




V 

D 




□ 


8 




A 


□ 




□ 


□ 


A 
□ 

A 


□ 




□ 


□ 




A 



v/(v ,Z ) 

^ nl p ' 



10 



FIG. 2: Ionization cross sections of hydrogen by fully stripped ions. The scaled experimental data 
are from Fig.l. Note that the data do not merge into a single curve. 

Resulting plots of the scaled cross sections are shown in FiglHl Comparing Fig.2 and Fig. 3 
one can clearly see that all of the experimental data merge close to each other on the scaled 
plot based on Eq. (j^ . 

The resulting universal function can be fitted with various functions, but the simplest fit 
was proposed by Rost and Pattard They showed that if both the cross section and the 
projectile velocity are normalized to the values of cross section and projectile velocity at the 
cross section maximum, then the scaled cross section a/a^nax is well described by the fitting 
function 

^/ ,N ^ exp(-t;Lx/^^ + l) /<.7^ 

a{V) = CTmax 771-^ • (.^7 j 



^max 



Here, (Jmax is the maximum of the cross section, which occurs at velocity fmax- For the 
present study (the case of the ionization cross section by the bare projectile), we predict 
that 

+ 1) III 



CTmax = TXalBni ^ 1) J^' ^^^^ 



t-max = f„zV^p + l, (29) 

where the coefficients B^i depend only weakly on the projectile charge. From Fig.3 one can 
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FIG. 3: Ionization cross sections of hydrogen by fully stripped ions showing the scaled experimental 
data and the theoretical fits. BA denotes the Born approximation [Ea. H21|) ]. Gillespie denotes 
Gillespie's fit according to Ea. ((^ . R.&P. symbolizes the fit proposed by Rost and Pattard [5^ in 
Eq.^. "New" denotes the new fit given by Eq.(|Hni)- 

estimate Bni = 0.8 for the ionization of hydrogen by protons, while for ionization of hydrogen 
by bare nuclei of helium and lithium, we find Bni = 0.93. As can be seen from Fig.3, the 
function in Eq.(j2Zj) with dmax and fmax defined in Eq.(|^Hj) describes well the cross sections 
at small and intermediate energies, but underestimates the cross section at high energies. 
The reason is that the function in Eq. ()27p does not reproduce the logarithmic term in the 
Bethe formula in Eq. ()ll|) . To improve the agreement with the experimental data and the 
Bethe formula we propose a new scaling for the fitting function in Eq. ()2f)|l defined by 

Qne.n^^^ ^ ^^P(-Va^') ^^ ^6 + 0.283 lu (2x2 + 25)] . (30) 

At large x >> 1, Eq. (j30|) approaches the Bethe formula in Eq.(|15|). and at small x < 1, 
Eq. ()30j) approaches the result in Eq. ()27j) . The function G"''^"'(x) has a maximum at x ~ 1, 
with ~ 0.86. Because 0.86 is in between the maxima of the scaled cross section of 

hydrogen by protons {Bni = 0.8) and the cross section for ionization of hydrogen by bare 
nuclei of helium and lithium {Bni = 0.93), we did not incorporate the coefficients Bni in 
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Eq.()3Up. This gives it a general form and introduces small errors of less than 8%. 

We have applied the new fit in Eqs. ()26p and (jHUj) to the ionization cross section of helium, 
shown in Fig. 4a. The symbols in Fig.4a denote the experimental data for if"*", He~^^, Li^^ 



53, 



5J, for C+6 [55|^ for I+^v and U+^p |56|, and for Au+^p [57], where = 10 - 40. The 



solid curves correspond to the continuum-distorted-wave-eikonal initial state (CDW-EIS) 
theoretical calculation from Ref. {5^, which is a generalization of the Born approximation. 
The CDW-EIS theory accounts for the distortion of the electron wave function by the pro- 
jectile. From Fig.4a it is evident that the CDW-EIS theory overestimates the cross section 
near the maximum, and underestimates the cross section at small energies. 

Direct application of the scaling formula in Eq.(|25p to the ionization of helium does not 
produce similar good results to the hydrogen case [see Fig. 4(b)]. But after applying the 
new scaling in Eq. ()26p . all of the experimental and theoretical results merge close together 
on the scaled plot, as is clearly evident in Fig. 4(c). Moreover, if we use the fit function 
of velocity normalized to the orbital velocity Vni estimated from the ionization potential of 
helium {I He = 24.6eV) making use of Eq.Q, the cross section is given by the same scaling as 
in Eq. ()2fi|l with the same function as in Eq. ()30|l . as evident from Fig. 4(d). (The number of 
electrons in the helium atom is Nni = 2, and therefore the scaled cross section is twice that of 
hydrogen.) From Fig. 4(d) it is clear that the new proposed fit in Eq.()26|) using the function 
in Eq. ()30p gives very good results for both hydrogen and helium. Further verification of the 
new scaling is difficult because reliable experimental data and numerical simulations for a 
broad range of projectile velocities are absent for other atoms. The discrepancy between the 
new fit and the helium data at very small velocities is discussed in the next section. 

Note that one experimental point in Fig.4 for C~^^ projectiles is located far away from 
the fit. The error bar for this point is about 30% 55|. This data may be inaccurate, as the 
experimental point is higher than the predictions of CDW-EIS theory, which overestimates 
the cross section near the maxima of the cross sections for all other ions. The reason for the 
large scatter in the uranium data on the scaled plot at small energies is not clear, because 
the experimental data for all other projectiles are located much closer to the fit line. 
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Exp □ He'' V Li'' A c'' Au'^' + f^^ U'^' 
CDW-EIS H- - He'' • • Li" O"; Fit new 

(a) 




E/(Mp(Zp+1 )) (keV/amu) EI{M^{Z+^) (keV/amu) 

FIG. 4: Ionization cross sections of helium by various stripped ions. The sohd curves correspond 
to the CDW-EIS theoretical calculation, and the symbols label the experimental data (see text for 
details). Shown in the figures are: (a) the raw data; (b) the scaled data from Fig. 4a, making use 
of Eq. (|25|) : (c) the scaled data making use of Eq. (|26|) : and (d) the experimental data scaled using 
only Ea. H26() together with the fit function. The notation "new fit" denotes Eg.ipU)). 

IV. THEORETICAL JUSTIFICATION FOR THE NEW FIT FORMULA FOR 
IONIZATION CROSS SECTION 

In this section we discuss the theoretical foundations for the new fit to the ionization 
cross section given by Eq. ()2(i|l and Eg. We start with an analysis of high projectile 
velocities. 
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A. Behavior of cross sections at large projectile velocities v > Vni 



In the region of high projectile velocities the new fit predicts the ionization cross section 



7)4 Z2 

high-energy, _ . 2 "^0 P 



0.566 In , + 1.26 



(31) 



which differs from the Bethe formula in Eq. pHl . [The factor a/ {Zp + l)/2 appears in the 
denominator under the logarithm in the first term on the right hand side of Eq. (jHT|) .] We 
claim that incorporating this factor gives a better cross section estimate than the Bethe 
formula. A comparison of the existing experimental data with the Bethe formula in Eq. ()ll|) 
and the fit formula in Eq. ()3ip is shown in Fig|3 The experimentally estimated uncertainty 
of 5.5% 42] is shown by the error bar. The region of validity of the Born approximation 
and, hence, the Bethe formula is 

V > max{2ZpVo,Vni)- (32) 



The first condition in Eq. ()32|) assures that the projectile potential is taken into account 
in the Born approximation; the second condition allows use of the unperturbed atomic 
wave function. Unfortunately the experimental data exists in the region in Eq. ()32|l only for 
the ionization of hydrogen by protons. Figure El shows that the Bethe formula describes 
the experimental data for ionization of hydrogen by protons within the error bar only for 
V > Gvq. Application of the fit formula instead of the Bethe formula reduces discrepancy 
with the data. 

The applicability of the Born theory and the Bethe formula in Eq. pij) was studied 
experimentally in Refs. It was confirmed that the necessary condition for 

the validity of the Bethe formula is given by the condition in Eq. ()32|) . The failure of the 
Bethe formula for large Zp is apparent from the experimental data for gold ions shown in 
Fig. 4(a). The ion velocity corresponds to v = 12vq ot v = 8.9vni, whereas Zp = 24,43,54, 
and does not satisfy the condition in Eq. ()32|) . As a result, the cross sections are much smaller 
than given by the Bethe formula, as evident from Fig. 4(a). (At large projectile energies, all 
data merge to the Bethe formula, which corresponds to a straight line in a logarithmic plot, 
similar to Fig.l.) 



The applicability of the Bethe formula is limited by the validity o: 



' the Born approxi- 



mation. One of the easiest ways to correct it was suggested in Ref.|6J|. Firstly, the Born 
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FIG. 5: Ratio of ionization cross sections of hydrogen by fully stripped ions to the Bethe formula in 
Eg. pij) and the fit formula in Ea. H31j) at high velocities. The experimentally estimated uncertainty 
of 5.5% is shown by the error bar. 

approximation is considered, making use of a classical trajectory for the projectile and a 
quantum mechanical description in the Born approximation for the electron. In this approx- 
imation, the probability of ionization or excitation is a function of the impact parameter p. 
Here, for brevity, we shall consider only ionization of the hydrogen atom. The projectile par- 
ticle interacts with the atomic electron with a potential energy V(R, r^) = — Zpe^/|R — re|, 
where R(t) = p + vt is the classical trajectory of the projectile particle, and re describes 
the position of the electron relative to the nucleus of the atom. For any impact parameter 
p, the probability of ionization is given by the square of the transition amplitude 

(33) 

Here, AE is the transferred energy in the transition, and \E'j and \E'/ are the initial and 
final electron wave functions, respectively. It can be shown that the calculations of ion- 
atom ionization cross sections using the conventional Born approximation describing the 
collision making use of momentum transfer (outlined in Appendix B) and the semiclassical 
Born approximation making use of the assumption of the straight line classical projectile 
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Pba{p) 



1 



trajectory [Eq. (jHH|) ] are equivalent 

For large impact parameters p >> ao, we can expand \^(R, r^) in powers of Veni/R 
according to 

nR,rJ = V^(-^ + ^). (34) 

The first term does not contribute to the matrix element in Eq. ()3Hj) due to the orthogonality 
of the final and initial states. Substituting Eq.(|n^ into Eq.(j2ni) and integrating in time 
yields 



Pba{p) 



2ZpUo 
pv 



rfre^i(re)^}(re) 



up J. (ujp 

XeKi 

V \ V 



+ IZe Kq 

V \ V / 



(35) 



where u = AE/h, and Kn is the modified Bessel function. Expanding the Bessel functions 
for small and large arguments, or simply evaluating the integrand in Eq. approximately, 
we can approximate 

^ 1, ^ <1 



0, 



(36) 



and neglect the second term on the right hand side in Eq. (jH3j) . which is small compared 
with the first term. The probability of ionization vanishes for p > pmax — v/uj = 2aof/fo, 
corresponding to the adiabatic limit. For p > Pmax, the collision time Pmax/'^^ > clq/vq is much 
longer than the electron circulation time around the nucleus, and the collision is adiabatic. 
Consequently, the ionization probability is exponentially small for p > 2aQv/vQ. 

The square of electron dipole matrix element averaged over all possible momenta of the 
ionized electron is 13811 



rJl^ = 0.283a^. 



(37) 



Note that the sum over all final states including both ionization and excitation gives 

^ < 0|se|/ >< /|se|0 >=< 0|x^|0 >= ^ < 0|r,2|0 >= al. (38) 



In this sum, 0.717 corresponds to excitation, and 0.283 corresponds to ionization 

For large impact parameters the momentum transfer to the electron is small and we 
can neglect the electron kinetic energy of the ejected electron compared with the ionization 
potential. As a result, AE ^ Ih = Eq/2 and uj = vo/2ao (in atomic units). Finally for 
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p > ao, the ionization probability is 



« 0.283 




1, p < 2aov/vo 
0, p > 2aov/vo 



(39) 



The ionization cross section is given by the integral 



poo 



a = 2n 



/ PBA{p)pdp. 



(40) 



For p > ao, we can use Eq. flS^ to estimate Pba{p)- For p < a^, the dipole approximation in 
Eq.(jnH) is not valid. To evaluate Pba{p) approximately for p < ao, we can utilize the fact 
that / (ite*^^*/'^V(R, r^) is a weak function of p for p < ao, and therefore Pba{p) ~ PBA{cio)- 
Substituting Pba{p) ~ -Pba('^o) ioi p < Qq, and Pba{p) from Eq. ljH^ for p > ao, into Eq.(jin)) 
gives 



The first term in Eq. ()41|) comes from contributions of impact parameters p < ao, and the 
second term originates from contributions of large impact parameters p > ao, respectively. 
Comparison with the exact result in the Born approximation in Eq. ljlH) shows that the 
contribution of impact parameters p < ao is underestimated, and 1/2 should be replaced by 
1.52. The above considerations are valid if the total probability of ionization and excitation 
[Pba{p) = {2ZpaoVo/ pv)'^ , for p > ao] for the entire region of impact parameters is less 
than unity, which requires 2ZpVQ/v < 1. (Note that the total probability of ionization and 
excitation is about 4 times larger for ionization only.) 

For 2ZpVQ/v > 1, the total probability of the ionization and excitation P^Xip) calculated 
using the Born approximation is more than unity, PbaIp) > -'-) f*^^ impact parameters 
P < Pbreak = 2ZpaoVo/v, indicating the breakdown of the Born approximation j6^. Similar 
to the previous case, we can estimate the ionization probability Pba{p) from Eq. ()39|) for p > 
Pbreak > o-o cLud assumc Pba{p) ~ PBA^Phreak) = 0.283 for p < Phreak- Thcsc Considerations 
result in a cross section estimate similar to the Bethe formula but with the logarithmic term 
in the form ln(pmax/Pmin) = ln(f^/fo^p), which gives 



This calculation results in a smaller cross section than the Bethe formula for 2ZpVQ/v > 1. 
Note that in the above analysis we have used unperturbed electron wave functions, which is 
valid only for w >> t>o. 




(41) 




(42) 
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While a number of smart semi-empirical ways to improve the first Born approximation 



were developed 



6G 



61 



62], the rigorous approaches to improve the Bethe formula are 



based on the eikonal approximation instead of the Born approximation 63[. The eikonal 
approximation is justified if fca„/ > 1, where k is the projectile particle wave vector k = 
Mv/h, and the projectile kinetic energy is large compared to the potential energy interaction 
with the target. For heavy projectile particles with mass much larger than the electron mass, 
these conditions are well satisfied. The ionization cross section in the eikonal approximation 
is given by 



a 



2tx 



qdq 



l/(q)r, 



where /(q) is the amplitude of ionization with momentum transfer q 

'i J Vdz 



/(q) 



k 

2TTi 



pdp < final \ exp 



h 



iq ■ p \initial > . 



(43) 



(44) 



The eikonal approximation in Eqs. ()43p and (j44j) accounts approximately for all orders of 
the perturbation series, whereas the Born approximation only make use of the first order. 
The calculations in the eikonal approximation yield a formula similar to Eq. ()42p 6^. Note 
that the validity of the eikonal approximation in Eg . (j44p is limited to v » fo,because the 
electron wave functions \E'j and are assumed to be unperturbed atomic functions. The 
influence of the projectile on the electron wave functions has to be taken into account for 
V < Vq. This is typically performed in the distorted wave approximation [16]. 

Therefore, the correction to the Born approximation in Eq. ljl^ and the eikonal approxi- 
mation give a formula similar to Eq. (j3T| but with a factor a^fT^p (a is a coefficient of order 
unity), instead of a/ (^p + l)/2. At large velocities, both formulas give similar results. 



B. Behavior of cross sections at small projectile velocities v < Vni 

If the projectile velocity is small compared with the orbital velocity, the collision is 
adiabatic and the electron circulates many times around both nuclei. The electronic energy 
states need to be determined in such a quasimolecule as a function of the positions of both 
nuclei at a particular time. In both the quantum mechanical and the classical approaches, 
ionization is only possible if during the collision the initial and final electronic terms cross at 
some instant. In classical mechanics this corresponds to the so-called "f/2 mechanism". In 
a coUisional system comprised of two nuclei of equal charges (say ionization of hydrogen by 
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a proton), an electron which is exactly in between the two nuclei experiences a very small 
electric field because the electric fields from both nuclei exactly cancel for all times at this 
point. The electron can "ride" this saddle point of the potential if its velocity is equal to 
one-half the velocity of the projectile. The collision dynamics is illustrated in Fig. 6. 



trajectory during ionizing collision 





FIG. 6: The trajectory of a v/2 colUsion is shown in Figs. 6(a) and (b). The initial conditions 
correspond to a hydrogen atom with total energy —1/2, and at f = —60 x = = y,Vx = = Vy, 
z = —1.606756 (solid line) and z = —1.606751 (dotted line). The projectile moves along z = 1 
with velocity 1/2. Atomic units are used: velocity is normalized to vq; distance is normalized to 
ao; and time is normalized to ao/fo- Figure 6(b) shows the position [x{t), z{t)] of the electron as a 
function of time, and the distance between the electron and the first (pi) and the second proton 
(/92) for the same conditions as in Fig.6(a). The trajectory of a S'-promotion is shown in Figs. 6(c) 
for fixed positive charges (v ^ 0). The initial conditions correspond to an internuclcar separation 
2ao (in atomic units), initial position of the electron z = 0, x = 1; and initial velocity Vx = 0, 
Vz = 1.155 (solid line), and Vz = 1.165 (dotted line). 
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From Fig. 6 one can see that the electron is stranded in between the protons at t = 15ao/fo 
and its velocity projection on the x-axis is one-half of the projectile velocity. A small 
variation of the initial condition from z = — 1.606756ao (solid line) to ^ = — 1.606751ao 
(dotted line) completely changes the result of the collision. After the collision the electron 
stays near the first nucleus and does not become ionized. As a result, the probability 
of ionization is extremely small even though the projectile velocity is not small (for the 
conditions in Fig. 6, v = 1/2 in atomic units). The mechanism for ionization described 
above is also so-called T-promotion in quantum mechanical descriptions j67| . 

Another mechanism for ionization is attributed to the so-called S-promotion mechanism 
[g?]. It is associated with the special type of trajectory of the electron in the field of two 
positive charges, shown in Fig. 6(c). Figure 6(c) shows that an electron with particular initial 
conditions tends to spiral with a large number of turns enclosing a segment of the straight 



line joining the nuclei Fig. 6(c) [68|. Such a trajectory is unstable - a small variation of initial 
conditions results in a completely different trajectory as shown in Fig.6(c). Analysis of the 
electron motion in the field of two positive charges, Zt and Zp, which are separated by a 
distance R is best described in elliptical coordinates 

where and tt are the distances from the electron to the projectile and target nuclei, 
respectively. Making use of atomic units, the classical trajectory in terms of the variables ^ 
and 1] can be expressed as 

_ 4(^^ - l)Pg dr] _ A{r]^ - 1)P^ 
where the canonical momentums and are 



(46) 



1/2 



.,^0,.,.^---f^-^)'^ 

Here E < is the total energy of the electron, P^ = C,Tid(f) / dt is the rotational momentum 
around the straight line joining the nuclei, and A is the integral of motion (for stationary 
nuclei) 

X = - — (P^ + + R{Zp cos 9p + Zt cos 9t). (49) 
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Here, ( is the closest distance from the electron to the straight line joining the nuclei; 



total rotational momentum; and Bp and 6t are the angles between Vp and R, and and 
— R, respectively. Moreover, Fp is the radius vector from the projectile to the electron; r^- is 
the radius vector from the target nucleus to the electron; and R is the radius vector from the 
projectile to the target nucleus. The canonical momentum in Eq. ljTTj) tends to infinity if 
,^ ^ 1, preventing the electron from approaching a segment of the straight line joining the 
nuclei, ^ = 1. In the special case 



the singularity vanishes at the point ,^ = 1 in Eq. ()47j) . As a result, for initial conditions 
satisfying the condition in Eq. ()5()|l . is finite for ^ = 1. From Eq. ()4fi|l . ^ approach unity 
exponentially with time - the limiting electron trajectory lies on the internuclear axis - as 
shown in Fig.6(c), where the initial conditions for the solid line correspond to the condition 
in Eq. ()50|) . A small departure from the condition in Eq. ()50|) shown by the dotted line in 
Fig. 6(c) prevents the trajectory from approaching ^ = 1. Thus the internuclear axis ^ = 1, 
represents the locus of points of unstable equilibria. In a quantum mechanical treatment, 
such periodic unstable trajectories is responsible for S-promotion of electron to the the 
continuum (ionization) when the nuclei approach each other [^. The potential barrier in 
Eq . ()47|) increases when R decreases. As a result, an electron near the top of the barrier slows 
down and is then collected and promoted to the continuum as the top of the barrier further 
rises. Due to the strong instability of the locus, a numerical simulation of the corresponding 
classical trajectory is extremely difficult. [We could not present the classical analog of the 
ionization scenario for S-promotion, in contrast to the T-pro motion as shown in Fig. 6(a) 
and (b).] 

The probability of ionization is greatly enhanced in quantum mechanics due to tunnelling 
into classically forbidden regions of phase space. The cross sections can be calculated using 
the quasiclassical method, where the probability of transition is given by 



is the vector dot product of the electron momentum with the C-axis; 



(r X p)^ is the 




(50) 




(51) 



where 




(52) 
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Here, S{p, e) is the classical action of the projectile ion, and p = \/2M[e — U{R, p) — Ei) is 
the projectile momentum, generalized to classically forbidden regions of phase space where 
p is complex 0]. The integration contour in Eq. (j5ip is in the complex R plane around the 
branch points {R'n) where the initial and final electronic terms cross [Ef{R^ = Ei{R'^)]. 
Moreover, n numerates different branch points or channels of ionization for S and T- 
promotions. The resulting cross section for hydrogen ionization by collision with a proton 
is Ifill 

0-adiabatic{v) = TTW ^ R^e''^^''^'" , (53) 
n 

where n labels many different channels, and the coefficients A„ and Rn are of order unity 
in atomic units {Rn is determined by the branch points -Ren)- In the range of projectile 
velocities v = 0.4 — 1, we find that Eq. (jHH|) can be approximated to within 10% accuracy 
by only two exponents with Ri = 1.9 , Ai = 0.53 (corresponding to S-promotion) and 
i?2 = 6.7, A2 = 1.8 (corresponding to T-promotion) . Because Ai << A2, primarily the 
S-promotion determines the ionization cross section at small velocities {v < 0.5), while both 
mechanisms contribute to ionization for v in the range f = 0.5 — 1. Recent experimental 
study and quantum mechanical calculations using the continuum-distorted-wave eikonal- 
initial-state (CDW-EIS) model jl^ show that a electron emission spectrum is dominated 
by a well defined electron capture to continuum (S-promotion) peak although existence of 
saddle-point electron emission (T-promotion) is not confirmed. 

The new fit predicts an extremely small cross section at very low velocity crj°^~^"^'"^^(f ) ~ 
exp(— l/t>^), whereas Eq. (j3^ gives cr adiabatic{v) ~ e~^'°/''. Therefore, the numerical fit in 
Eq.(jnni) underestimates the cross section for v < 0.5, but gives a result close to the sum in 
Eq. ()53j) for v in the range v = 0.5 — 1. While the data for hydrogen at very low projectile 
velocity is absent, and the fit agrees well for the entire dataset in Fig. 3, the disagreement 
is clearly seen when the fit is compared with the experimental data for the ionization of He 
shown in Fig. 4(d). Adiabatic theory results are absent for helium, but the experimental ion- 
ization cross section of He by protons can be described by Eq. (|53p with different coefficients 
A„ and Rn- The behavior of the experimental ionization cross section of He by He"*"^ is 
somewhat puzzling because of the very slow decrease of the cross section for small projectile 
velocity. 

In view of these observations, the applicability of the new fit is limited to 
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V l\oni\/ {'Zp + 1)] > 0.5. Note that for small projectile velocity the ionization cross sec- 
tion is ten times smaller than the maximum of the cross section, cxmax, and the ionization 
cross section is completely dominated by charge exchange, whose cross section is compara- 
ble to (Tinax- Consequently both experimental measurements and theoretical simulations are 
very difficult for very small projectile velocity. 




V. CONCLUSIONS 

The new scaling in Eq. ()26|) for the ionization and stripping cross sections of atoms and ions 
by fully stripped projectiles has been proposed. The new scaling does not have any fitting 
parameters and describes the shape of the cross section as a single function of the scaled 
projectile velocity [Eq. (p?njl ]. Note that previous scaling laws either used fitting parameters 
1) or actually did not match experiments in a wide range of projectile velocities 
The proposed scaling formula agrees well with theoretical predictions in the limit 
of large projectile velocities. The new scaling has been verified by comparison with available 
experimental data and theoretical simulations for the ionization cross sections of hydrogen 
and helium by i/'*', i/e"'"^, Li"*''^, and O^^ . The agreement between the new proposed 
scaling and experimental data is very good. The difference between the proposed fit and 
the experimental data is within 15% accuracy, which is similar to the estimated uncertainty 
in the measurements. The validity of the fit is limited at very small velocities, where the 
ionization cross section is very small, about one-tenth of the maximum cross section o"max, 
and the ionization cross section is completely dominated by charge exchange, whose cross 
section is comparable to cXmax- Finally, the fit is valid for scaled projectile velocity v > 



0.5vni\/Zp + 1, where Vni = vq a/ 2/„/ / Eq is the orbital velocity of the electron estimated from 
the ionization potential where Eq = 27.2eV (twice the hydrogen ionization potential). 
Similarly, the fit is valid for E > 12.5{Zp + Vjlni/E^ in units of keV/amu, where E is the 
projectile kinetic energy per nucleon. 
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APPENDIX A: CLASSICAL CROSS SECTION AVERAGED OVER ATOMIC 
ELECTRON VELOCITY DIRECTIONS 



Gerjuoy averaged the Rutherford cross section over all orientations of the electron velocity 
Ve (for a fixed electron speed Ve) and derived the differential cross section da /dAEjvp, y, AE) 



for energy transfer AE in the collision between a free electron and the projectile j3J]. The 
total cross section is calculated by integrating over values of energy transfer larger than 
the ionization potential {AE > Ini ) and averaging over the electron velocity distribution 
function (EVDF) /(fe). This gives 



(T{v,Inl,Zp) = Zp / ai^X'^^^e) f {Ve) dVe 



where 

and da/dAE{ve,v, AE) is defined by 



da 
dAE 



[v,Ve,AE)dAE, 



(Al) 



(A2) 



^, A P\ ™o ^0 Siv, Ve, AE) 



dAE 



4 AE^ v'^Ve 



(A3) 



where 



S{v,Ve,AE) 



{v' - vl) {vl -v^~ 2AE/me) {v^^l - <) + 
2 {vl + v^ + AE/rrie) {v^^ - vi^^) - 1/3 {vl^ - vlJ) 



Here, Vup and viow are defined by 

Vup =Ve + V, 

Vlow = max (^\ve - v\ , ^/ v"^ - 2AE/me -v]. 

For very large projectile velocities v >> v^, it follows that S ~ 8ve (2fg/3 + AE/nie), and 
Eq.(jX3} yields 



da 



high energy 
classical 



{v, Ve, AE) = 27ra, 



2 ^0 



2mev! ^ 



(A4) 



dAE ' " °AE3met;2 V 3 

Substitution of Eq. ()A4|) into Eq. ()A2|) , and subsequent substitution of Eq. ()A2|) and the EVDF 
Eq.® into Eq-dH} give 



high energy/ j \ ry2 2 '^O'^U 

^classical l"^) -'ni, ^pj - ^TTZ. ag-— — . 
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vlE^ 



(A5) 
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In the general case with v ~ We, substituting the EVDF Eq.® into Eqs. ()A2|) and ()A1|1 yields 



^■2 

(^classicali^'V^ Ini, Zp) = TTCLqEq—^G classical ( / = ) ; (-^6) 



nl 

where 



Gclassicalix) = ^ / T"^;:^ dAEdVe. 

Jo Jl/2 AE^Ve 



The approximate formula for G ciassicai{x) is given below in Eq. ()C3|) . 



APPENDIX B: THE BORN APPROXIMATION 



Although the Born approximation is valid only for large projectile velocities v » ZpVo 
|27|. the Born approximation does give results close to the experimental data even outside 



its validity range 



3- 



Therefore, we have studied cross sections in the Born approximation 



for the entire velocity range. 

In the Born approximation, the ionization cross section for 
fully stripped projectile atoms with charge Zp is given by 



lydrogen atoms by impact of 



38, 



aZ\v) = Malzfl r ^^-^dq, (Bl) 

Jo H 

where P/^, (g, v) is the probability of ionization, and qrrieVQ is the momentum transfer dur- 
ing the collision. We introducing the velocity in atomic units v = v/vq, and Pj^^{q, v)is 
determined by js^ 

Here, 0(x) is the Heaviside function, and dP{q, K,)/dK is the differential probability of eject- 
ing an electron with momentum KnieVQ when the momentum transfer from the projectile is 

qrrieVo, 

= \mp)MP + ^))\" = \^^^^^ (B3) 
In Eq. ()B3|) . ^I/*(p) and ^^^^(r) are the wave functions of the continuous spectrum (ionized 
electron) in momentum space and coordinate space, respectively; ^I^o(p) and \I^o(r) are the 
wave functions of the ground state, and star (*) denotes complex conjugate. According to 

dPjq, k) ^ ^8^ , [q^ + 1(1 + K^)] exp{-2//tarctan[2/t/(l + q' - k^)]} 

dn ''^ [(g + + 1]3 [^^ _ + 1]3 (1 _ 
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For q » 1, the function dP{q, K)/dK has a sharp maximum at k = g j2 

dP{q, k) 8 1 



(B5) 



37r[(g-fc)2 + l]3' 

which simply means that the entire momentum q is transferred to the ionized electron 
momentum k. At small q < 1, dP{q, K,)/dK, ^ nq^ and the width of the function P(g, k) as 
a function of k is of order unity in atomic units. 

For large projectile velocity v » vq, considerable simplification can be made by neglect- 
ing the electron kinetic energy in the argument of the Heaviside function in Eq. ()B2|) . 
The approximation 

e -e g-f (B6) 

V VO J \ VO J 

is referred to as the close-coupling approximation. In this case, P{q,v) can be characterized 
by a function of one argument, Smhig), with 

PiJq,v) = SUQ)Q{q-^-^y (B7) 

where 

S^nH{q) = r ^^^^dn. (B8) 
Jo dn 

The function Sinhio) is refereed to as the total ionization transition strength j46|. Sub- 
stituting Eg. ()B6|1 results in artificial, additional contributions to the integral in Eq. ()B2|l 
for K > Kadd = A/2(gi;/fo — Ini/ Eq). For large projectile velocities v » Vq and g >> 1, 
i^add — ^/'^Qv/vq. The function dP{q, k) / dn has a sharp maximum aX k = q [see Eq. ()B5|) ]. 
Therefore the artificial additions for k > Hadd do not contribute to the integral if Hadd > Q, 
which corresponds to g < 2t>, and the substitution in Eq. ()B6|) is valid. In the opposite case 
of large projectile velocities v » Vq but small g, it follows that g ~ volni/ (vEq) « 1 , for 
the range of g Hadd ~ 1, and the function dP{q, K)/dn decreases rapidly for k > 1. Therefore, 
the artificial additions for n > Kadd do not contribute to the integral if Kadd > 1- Hence, the 
substitution in Eq. ()B6|) is vahd for t> >> vq. Figure 7 shows plots of Pi^^{q,v) [Eq. ()B2|) ] and 
Sinhil) [Eq. ()B8|) ] for = 1 and v = ?>. At small projectile velocities v < fo, the substitution 
in Eq. ()B6j) produces a considerable error [see Fig. 7]. For repetitive calculations, the function 
Sinhil) in Eq. ()B8|) can be approximated to within 3% accuracy by 



(g-0.9)2+1.21 y ^ ^ 

tanh(0.8g) g > 2 



(B9) 
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FIG. 7: Total ionization transition strength for atomic hydrogen as a function of transferred 
momentum q. The exact function P{q,v) [Eq.( IB2|) ] for v = 1 and -u = 3 is compared with the 
approximate function Sinhiq) [Eq. (|B8|) ] (which is independent of v) and the fit S^J^{q) in Eq. (|B9|) . 

The functions Sinh{q) [Eq.dHH)] and Sf^Hq) [Eq.(|B9}] are shown in Fig. [7| 

Having estimated the function Pj^^{q,v), the total cross section can be evaluated ana- 
lytically for large v >> Vq. The region of small q contributes significantly to the cross 
section [see Eg. (jBll) ]. Therefore, we split the integration in Eq. ljBlll into the two re- 
gions q < qup and q > qup, where qup = 1/2. In the first region q < q^p, it follows that 
Pinii^j'^) ~ ^inhil) ~ 0.283g^, and the integration in Eq. ()Bl|) gives 

/ dq ^ ^ / dq- = 0.283 ln(g„p/g^i,), (BIO) 

Jo Q Jq^in ^ 

where gmin = VoIni/vEo. In the second region, only the range of qup < q < 2 contributes 
to the integral, because at large q » 1, Pi^^{q,v)/q^ ~ 1/q^ and the contribution to 
the integral for large q quickly decreases to zero. At very large q > 2v, Pi^^{q,v) became 
smaller than unity, but this region does not contribute to the integral and can be neglected. 
As a result, the integral dqPf^^(q,v) /q'^ does not depend on v (for the large v under 
consideration). The integration from q^p to infinity gives dqPi^^{q,v)/q^ ~ 0.666, and 
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finally the result is similar to the Bethe formula in Eq.()ll|) with 



^2 

^Bethe^-^ = Snal ^ [0.283 In (v) + 0.666] 



(Bll) 



The small differences from the Bethe formula are due to utilization of the close coupled 
approximation in Eg. (jBlOj) . which overestimates P7^,(g,f) at small q [see FigEj. 

Comparison with the exact calculation (Fig.l) shows that the Bethe asymptotic result is 
close to the exact calculation in Ea. llBlll for > 2. To extend the Bethe formula to lower 
velocities, the second-order correction in the parameter Vq/v has been calculated in |44| . 
yielding the cross section in the form 



Bethel 
mod \ 



Alia, 



~2 



0.571n (?^) + 1.26 - 0.66. 



(B12) 



where v = v/vq. Equation ()B12|) agrees with the exact calculation in Eq. ()Bl|) to within 
10% for V > 1.1. We have developed the following fit for the cross section in the Born 
approximation, 



47ial ^ [0.283 In (^J^ + l) + 1.26] exp 



1.95 



(B13) 



which agrees with the exact calculation in Eq. ()Bl|) to within 2% for v > 1, and to within 
20% for 0.2 <v< 1. 

The previous analysis was performed for the hydrogen atom. In the case of hydrogen-like 
electron orbitals, the similarity principle can be used. The quantity dP{q, n)/dK, is identical 
for different electron orbitals if q, k are scaled with the factor 1/ Zt = v^/vni [27]. Therefore, 
Pniilyv) = PniQVo/vnhv/vni), whcrc H denotes hydrogen atom, and 



V 
Vnl 



Ana 



oA 



[0.283 In (?^2 + l) + 1.26] exp 



1.95 



1^(1 + 1.2^^2 



where 



V 

Vnl 



(B14) 



(B15) 



^J21nllme 

As we have noted for helium, most scalings can be used even for non-hydrogen-like electron 
orbitals, provided the relationship in Eq. ()B15|) is used. 



1. Comparison between the quantum mechanical and classical trajectory calcula- 
tions for V » Vnl 

We have previously noted that the classical trajectory calculation underestimates the 
ionization cross section at large velocities v >> Vni- To compare the ionization cross section 
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q (a.u.) 




q (a.u.) 



FIG. 8: Probability of ionization of atomic hydrogen as a function of transferred momentum; Pc{q) 
is given by classical mechanics [Ea. ()B30|) ]. and Pq{q,v) is given by quantum mechanics [Ea.f IB2|) ]. 
The plots correspond to (a) v = 5 and (b) v = 15. 



calculated in the classical trajectory and Born approximations, we present both cross sections 
in the form of Eq. ljBlj) . In the limit v » Vni, the momentum transferred to the electron 
during a collision with impact parameter p is given by Eq.((H), i.e., 



Qxip) = rUeAvxip) 



vp 



(B16) 
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where x— axis is chosen in the direction perpendicular to the projectile ion trajectory along 
the momentum transfer. Because v » Vni, the electron velocity is neg lected in Ea. (lBT6|l . 
In classical mechanics, ionization occurs if the energy transfer to the electron is more than 



nl- 



the ionization potential, [(rngVe + q) — m^v ]/2mf, > J, 

A small momentum transfer to the electron along the projectile trajectory qz{p) can be 
determined making use of the energy conservation. Due to conservation of the momentum, 
the momentum transferred from the projectile particle is —qz{p)- The projectile energy 
change is [(Mv — q)^ — M^f ^]/2Af = —vq^- Conservation of energy gives 

vq, = ^\[m,w, + q)2 - m^z;^]. (B17) 

In the limit v >> Ve, it follows that << qx, and consequently the total transferred 
momentum to the electron is g = ^/q'^ + ql — qx- The momentum of the ejected electron 
can be determined from the energy conservation relation 

fi:V2me = [(meVe + q)' - mlvl]/2me - /„/. (B18) 

In classical mechanics, the ionization probability of the ejected electron with momentum 
K in a collision with total momentum transfer q is given by the integral over the electron 
distribution function, 

dPc{q,K) f nr ^ f <f r\ ^■R1o^ 

— / /(ve)dv/ — - - q^Vx -7r--Ini]. (B19) 



dn nie J ^ \2me 2m 

Introducing the one-dimensional electron distribution function 



Wex) = J f{ye)dVydVz, (B20) 



and substituting q ~ q^, Eq. ()B19|) simplifies to become 



dPc{q,K) _ 1^ J f /t^ - - 2mJni , ^^21) 



dn qrrie \ 2gme 

For hydrogen-like electron orbitals given by Eq.®, fx{vex) can be readily calculated to be 

fAVex) = ^ . ^ 5' 213 - ^^22) 

Kx + <i\ 

Substituting the hydrogen-like electron distribution function Eq. ()B22j) into Eq. ()B29|1 gives 
in atomic units 



dPc{q, k) IQk {2qme)^v, 



5/, ,5 
nl 



36 



3 • 



(B23) 



Let us compare Eg . (16231) with the quantum mechanical result Eq. ()B5|) . In the limit g >> 1, 
K, ^ q and the two functions are equivalent. Both functions dP{q, K)/dK have a maximum 
at K, = q, and the width of the maximum is of order 1, which simply means that the entire 
momentum q is transferred to the ionized electron momentum k. 

Moreover it is possible to prove that the classical mechanical dPc{q, K,)/dK is equivalent 
to the quantum mechanical function dPg{q, K,)/dK, for any s— electron orbital (spherically 
symmetrical wave function). Indeed, for large k » 1, the ejected electron can be described 
as a sum over plane waves ^I/*(r) ^ e^^^, and substituting ^I'*(r) into Eq. ()B3|) gives 

where integral over do^ = 27r sin i)di} designates averaging over all directions of the k- vector, 
is the angle between q and k, and / (ve) is the electron distribution function in velocity 
space. Note that |q — k|^= q^+k^ — 2q ■ k = {q — + 4qksmi)/2'^. In the limit q » 1, 
k q and only small d contribute to the integral in Eq. ()B24|) . Therefore, averaging over all 
directions of the k-vector gives 



1 /;fi^V^<'o..-^//KS±^|2rf«. (B25) 



Introducing v± = kd/me-, the integral in Eq. ()B25j) takes form 

' q — k 

nip 




(B26) 



where fx is the one-dimensional electron velocity distribution function. Substituting 
Eqs. frR26|l and (fR25|l into Eq.(lB24l) yields 

dPg{q, k) 1 j. f q — k 



, fA- ■ (B27) 

dK me \ nie J 

Note that in the limit q» meVni, it follows that n ^ q, and Eq. ()B21|) becomes 



dn nie \ f^e 

Finally, comparing Eqs. ()B27|) and ()B28|) we arrive at the equivalence of functions dP{q, k) / dn 
in quantum mechanics and classical mechanics in the limit q » nieVni- 

The situation is completely different for small q << nieVni ■ From Eq. ()B23|) it follows that 
dPc{q, n)/dK ~ nq^ , and dPc{q, K)/dn is much smaller than dPq{q, K,)/dK ~ nq"^. Therefore, 
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classical mechanics strongly underestimates the probability of ionization for small transferred 
momentum q < nieVni- 

The total probability of ionization in classical mechanics is 

P,(g) = r dn'^^^^i^ = [q (qv,, + /(ve)dv^. (B29) 







dK J \ 2m, 
Equation ()B29|) simplifies to become 



Pc{q) = y ® (^^^ex + - Ifil^ fx{Vex)dVex- (B30) 

The differential cross section for momentum transfer q is given by 

da^q) = 27rp{q)dp{q), (B31) 
where p(g) is given by Eg. ()B16)1 . Substituting p(g) from Eg. ()B16jl into Eg. ()B31)1 gives 

dadq) = -^T^dq, (B32) 

which is the Rutherford differential cross section for scattering at small angles. Finally, the 
total ionization cross section is 



a, = Svra^Zl^ j t^dq. (B33) 



'ini/v q^ 

In Eg. f)B33|) . we accounted for the fact that the minimum g is g = Ini/v. Note that 
in the region q = [1 — 3]Ini/v ionization occurs due the collisions with very fast electrons 
Ve ^ V » Vni, and qx ~ qz- The previous analysis which assumed Ve << v and qx » 
q^ is not valid in this region of extremely small q. However, because Pc{q)/q^ as 
g — > 0, this region of g = [1 — 3]Ini/v does not contribute to the integral in Eg. ()B33|) 
and can be neglected. Moreover such small momentum transfers correspond to very large 
impact parameter p/w ~ ani/vni, where the collision becomes adiabatic. Therefore, accurate 
calculations yield even smaller Pc{q) than in Eg. ()B30|) . 

Eguation ()B33|1 is identical to Eg. (jBlll . where the guantum mechanical ionization proba- 
bility Pq{q, v) is replaced by the classical mechanical ionization probability Pc{q) in Eg. ()B30jl . 
The functions Pg{q,v) [Eq.^] and P^{q) [Ea. (p30|i ; are shown in Fig. 8. Figure 8 shows 
that the functions Pi^^{q,v) and Pc{q) are nearly identical for g > 0.6. The classical proba- 
bility of ionization Pc{q) rapidly tends to zero for g < 0.6, while the guantum probability of 
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ionization, Pq{q) ~ 0.283g^, is much larger than Pc{q) at small q. The cross section is deter- 
mined by Pq{q)/q^- Therefore the region of small q contributes considerably to the quantum 
mechanical cross section. Note that Pq{q)/q^ ^ as g — Ini/v . It follows that the region 
of small q contributes most to the cross section [compare Fig.8(a) for v = 5, and Fig. 8(b) 
for V = 15]. For v = 5, the classical mechanical ionization cross section in atomic units is 
(Jc = 0.23, and the quantum mechanical ionization cross section is aq = 0.30, which is 30% 
larger than the classical mechanical cross section. For v = 15, Uc = 0.025 and aq = 0.043, 
which is 70% larger. 



APPENDIX C: FORMULARY FOR IONIZATION CROSS SECTION 

In the high energy limit of fast projectile motion v >> Vni, the classical mechanical 
calculation can be readily carried out (see Appendix A). 

The Bohr formula neglects the electron velocity in the atom completely, which 



gives 



a^"'^{v,I^i,Z,) = 27rZX^-^. (CI) 



Accounting for the electron velocity gives an additional factor of 5/3 compared with the 
Bohr formula. This gives the classical mechanical ionization cross section in the limit of 
high projectile velocity 

high energy/- j rj \ ^ 2 2^0-^0 

(^classical V" , ^nU Zp) = -2T:Zpa^y-—. 

"J V Inl 

In the general case with v ~ Vni, the classical mechanical calculation accounting for the 
finite electron velocity in the atom, but ne glec ting the influence of the target nucleus on the 
electron has been performed by Gerjuoy [2J] [see Appendix A]. This gives 

J-nl \ \/2Inllmr I 



The tabulation of the function G"-^*^^(x) is presented in Ref.jsoJ for x > 1, and in Ref . j37 1 
for X < 1, which gives 
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where 



f + e "ctanc + + 1. (35 - f - f ) + 

^ [(5 - 4x2) (3a2 + 1.5a6 + b^) - cx (7.5 + 9a + 56)] - 
i^xa^ ln(4x2 + 1) - ax^ (l + 2^e|si£) (2.5 + 3a + 4a2 + Sa^) 

a = l/{l + x^) c = 3x/4 6=1/(1 + 0^). 



(C4) 



and 



Gryzinski's approximation for the ionization cross section j20| expressed in the 
form of Eq. ()C2|) is given by 



a 



(C5) 



where 

G^'y^x) 



[a + |(l + ^)ln(2.7 + x)] (1-/3) + for a; > 0.206 



^x^ for X < 0.206. 



(C6) 



and a = x^jil + x^) /3 = l/[4x(l + x)\ 

Bethe's asymptotic quantum mechanical calculation in the Born approxima- 
tion [38] is valid for v/vq> 2Zp and v » Vni ^^"^ "^^^ t)e expressed as 



71472 

^Bethe ^O^P 



nl 







0.571n ( — ) 


+ 1.26 


\VnlJ 





(C7) 
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28] 



The region of validity of the Born approximation and, hence, the Bethe formula is 

V > max{2ZpVo,Vni)- (C8a) 

The first condition in Eq. ()(^8aj) assures that the projectile potential is taken into account in 
the Born approximation; the second condition allows use of the unperturbed atomic wave 
function. 

To describe the behavior of the cross section near the maximum, the second-order cor- 
.ection in the p^a^et. ha. been ca.cn.ated in Ref.Q. yie.d.ng the e.. seetion ,n 
the form 



Bethe 
^ mod 



?;2 



1 



0.566 In (?0 + 1.26 - 0.66— 



(C9) 
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where 



2 -^0 



Vnl \plT^m^ 

In the general case with v ^ Vni, the ionization cross section in the Born approximation was 
first calculated in Ref.j43]. We have developed the following fit for the Bates and GrifRng 
result 

= — 1 = 4vra^, 4 5 [0-283 In (^^^ + l) + 1.26] exp \- ^. 
The Bethe cross section valid for relativistic particles [39] is given by 

7n V„iZj„ 



(CIO) 



Bethe a^„2 "^0 ^nZ^p 



nl 2 



{ML,[21n(7A)-/5']+C^.on} 



(Clla) 



where [3^ = Vp/ c, c is the speed of light, 7p = 1/ ^/l — /5p, and Mf^^ and Cjon are characteristic 
constants depending on the ionized atom or ion. For the hydrogen atom, Mf^^ = 0.283 and 
Qon = 4.04. 

Gillespie's fit for the ionization cross sections j4a] is given by 

2" 



= exp 



-A„z (^t;oA/^/^) 



a, 



^Bethe 
mod 1 



(C12) 



where \ni is a characteristic constant of the ionized atom or ion (for example, for the ground 
state of atomic hydrogen, A„; = 0.76), and o'm'od^ modified Bethe cross section defined 

in Eq.jnni). 

The Olson scaling ^49^ for the total electron loss cross section which includes both 
the charge exchange cross section cr^*^ and the ionization cross section, is given by 



(j*''(f , Zj,) = na^ZpAnif 



Olson 



(C13) 



where /(x) describes the scaled cross sections 

1 



/ 



Olson I 



X) 



[l — exp (— x^)] 



and and Ani are constants. For example, 'Jh = a/5/4 =1.12 and Ah = 16/3 for atomic 
hydrogen, whereas = 1.44 and Ahe = 3.57 for helium. 

Rost and Pattard [52] proposed a fit for the ionization cross section, which utilizes two 
fitting parameters, namely the maximum value of the cross section and projectile energy 
corresponding to the maximum value of the cross section. They showed that if both the 
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cross section and the projectile velocity are normalized to the values of the cross section and 
the projectile velocity at the cross section maximum, then the scaled cross section a/cTmax 
is well described by the fitting function j52 1 

^(.A rr ^^P(-^max/^^ + l) /p. . N 

a(y) = (Jmax TT-^ , 1^14) 

V /V 
I max 

where (Tmax is the maximum cross section, which occurs at the velocity fmax- 

We have shown that for ionization by a bare projectile, the values (Xmax and fmax are well 
defined by the projectile charge Zp, with 

2 ^0 
^^max = Vnl\/Zj, + 1, (C16) 

where the coefficient B^i depends weakly on the projectile charge. For example, for ionization 
of hydrogen by protons, B^i = 0.8, and for ionization of hydrogen by bare nuclei of helium 
or lithium, Bni = 0.93. 

Equation ()(n4j) describes well the cross sections at small and intermediate energies, but 
underestimates the cross section at high energies, because it does not reproduce the logarith- 
mic term of the Bethe formula in Eq. ()C7|) . To improve the agreement with the experimental 
data and the Bethe formula, we propose the new scaling 

a-{v, Ini, Z,) = ^al-^§G-- ( ^==] , (C17) 

where 

Qnew^^^ ^ ^^P(-y^ ) [1.26 + 0.283 In {2x^ + 25)] . (C18) 

In all previous equations cross section are given per electron in the orbital. If Nni is the 
number of electrons in the orbital, the ionization cross section of any electron in the orbital 
should be increased by the factor Nni. 

Finally, it should be noted that a number of other semi-empirical models have been 
developed, which use up to ten fitting parameters to describe the ionization cross sections 
over the entire projectile energy range 
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